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Abstract 

In this work the charm and bottom quark masses are determined from 
QCD moment sum rules for the charmonium and upsilon systems. To 
illustrate the special character of these sum rules when applied to 
Coulomb systems we first set up and study the behaviour of the sum 
rules in quantum mechanics. In our analysis we include both the re- 
sults from nonrelativistic QCD and perturbation theory at next-next- 
to-leading order. The moments are evaluated at different values of 
which correspond to different relative influence among the theoretical 
contributions. 

In the numerical analysis we obtain the masses by choosing cen- 
tral values for all input parameters. The error is estimated from a 
variation of these parameters. First, the analysis is performed in 
the pole mass scheme. Second, we employ the potential-subtracted 
mass in intermediate steps of the calculation to then infer the quark 
masses in the MS-scheme. Our final results for the pole- and MS- 
masses are: = 1.75 ± 0.15 GeV, mdmc) = 1.19 ± 0.11 GeV, 
Mb = 4.98 ± 0.125 GeV and mb{mb) = 4.24 ± 0.10 GeV. 
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1 Introduction 



Quantum Chromodynamics, the fundamental theory of strong interactions, 
represents a basic building stone of the Standard Model. The determination 
of its parameters remains an essential task within modern particle physics. 
The strong coupling constant can be obtained from many different sources 
to rather high accuracy IJ. Investigations for the quark masses face much 
more severe problems. Confinement effects must be taken into account for 
most systems sensitive to the masses. Therefore, apart from the top mass, 
nonperturbative methods hke QCD sum rules lattice QCD (5j(6j or 

chiral perturbation theory [71IH1 have to be employed. 

The extraction of the heavy quark masses was among the first applications 
of the method of QCD sum rules [210] • In this framework the hadronic 
parameters can be related to a perturbative QCD calculation including the 
nonperturbative condensate contributions. The analyses were later updated 
and extended [Hllini- Recent times have seen a renewed interest in these 
investigations. With the development of nonrelativistic QCD (NRQCD) [TTl 
112] it has been recognised that the Coulombic form of the potential plays a 
major role in determinations of the charm and bottom quark masses. The 
contributions from NRQCD have been calculated up to next-next-to- leading 
order (NNLO) |13|i14li15. . In a region where the system is sensitive to mass 
effects, they dominate the theoretical evaluation of the sum rules over a pure 
perturbative expansion in the strong coupling constant. 

The fundamental quantity in this type of sum rule analysis is the vacuum 
polarisation function n(g^): 

n^.(g') =^Jd'x e^'^^ (T{>(x)jt(0)}) = {q,q. - g,.q') U{q') , (1) 

where the relevant vector current is represented either by the charm {x) = 
(c7^c)(x) or the bottom current j^{x) = (67^6) (x). Via the optical theorem, 
the experimental cross section a{e~^e" — >■ cc, bb) is related to the imaginary 
part of n(s): 

„ / N 1 cr(e+e^ ^ cc,bb) ^ „ , 

Usually, moments of the vacuum polarisation are defined by taking deriva- 
tives of the correlator at s = 0. However, in this work we will allow for an 
arbitrary evaluation point s = —Am'^^ to define the dimensionless moments: 

x„(f) = — (Wjj) nw 



(3) 
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As we will see in more detail later, the parameter ^ encodes much information 
about the system. By taking ^ larger the evaluation point moves further away 
from the threshold region. As already discussed in jHj, this leads to a better 
convergence for the perturbative expansion. The price to be paid is a small 
dependence of the moments on the mass. This again limits the possible 
accuracy when extracting the mass from the moments. Using a dispersion 
relation we can write the moments AiniO ^ integral over the spectral 
density R{s): 



V J J (s + 4m2^)"+i J (1 + ^(1 - 
The last equation represents a convenient way to express the moments through 



the velocity of the heavy quarks v = — Am? / s. 

The moments can either be calculated theoretically, including Coulomb 
resummation, perturbation theory and nonperturbative contributions, or be 
obtained from experiment. In this way one can relate the heavy quark masses 
to the hadronic properties of the quark-antiquark systems. The theoretical 
setup is identical for the charmonium and upsilon and in the first chapters we 
will thus not specify the quark content. The difference will become crucial for 
the phenomenological part and the numerical analysis. First we will perform 
the analysis for the bottom quark mass where the theoretical expansions 
converge better and then determine the charm quark mass from the more 
delicate charmonium system. 

Of decisive importance for the determination of the masses is the thresh- 
old behaviour. Here the system reacts very sensitive to mass effects. A small 
change in the mass leads to a relatively large variation of the moments. Thus 
the mass can in principle be determined with rather high accuracy. Therefore 
we try to develop a consistent physical description of the threshold region 
which includes all theoretical contributions and perform a matching between 
the low and high energy region. 

Before turning to a discussion of the individual contributions, the next 
section shall highlight the special character of the sum rules when applied to 
the Coulomb system. To this end, we present the sum rules in the framework 
of quantum mechanics which were developed in ^H] • Then we present here for 
the first time the application of the moment sum rules to a system governed 
by the Coulomb potential. All theoretical contributions can be described 
analytically. The sum rules show the dependence of the mass on the moments 
explicitly. Furthermore, the dependence of the pole and continuum parts on 
n and ^ will be studied and we investigate how these parameters determine 
the relative influence between the different contributions. 
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The results from NRQCD can be described similar to the quantum me- 
chanical sum rules. They contain a nonrelativistic Green's function which is 
composed of a continuous spectrum above and poles below threshold. We 
will isolate these contributions and analyse their influence on the mass and 
the error separately. Namely the poles will give the largest individual contri- 
bution to the mass. Therefore we will not use the nonrelativistic expansion to 
the energy levels itself which suffers from large corrections, but rather evalu- 
ate directly the Green's function in a region where the expansion is expected 
to work. 

As already done in our previous work for the charmonium system in 
the theoretical description of the correlator we will use information both from 
Coulomb resummation and perturbation theory. The perturbative expansion 
contains a part which is not included in NRQCD and will be added to the cor- 
relator. The nonrelativistic spectral density is valid for low velocities whereas 
the perturbative spectral density is well suited for high velocities. Therefore 
we will introduce a separation velocity to construct a spectral density for the 
full energy range. This will allow us to obtain a stable mass prediction for a 
wide range of parameters. The matching will turn out to be very smooth for 
the upsilon system but leading already to a gap for the charmonium. The 
analysis will show, however, that this mismatch has no numerical influence 
on the mass. The parameter ^ will be used to shift the analysis towards 
a more perturbative region. There the expansions of NRQCD and pertur- 
bation theory converge faster and the different theoretical contributions are 
more equally distributed. This will reduce the systematic uncertainty of the 
sum rule approach. Since the sum rules keep the analytic dependence on the 
input parameters, we will investigate in detail their influence on the mass. 
To estimate the error we will vary each of the parameters in suitably large 
chosen windows. The dominant uncertainty will come from the nonrelativis- 
tic expansion and manifests itself in a large dependence on the corresponding 
scales. 

Until now we have not specified the mass definition to be used in eq. (jH)). 
A natural choice for this mass is the pole mass M. In the numerical analysis 
we will first use the pole mass scheme to extract the pole masses. However, 
as the pole masses suffer from renormalon ambiguities ^H], they can only be 
determined up to corrections of order 0{Kqcd)- Therefore, in the second part 
of our analysis we shall use the potential-subtracted (PS-) mass mps [TH] . 
From this mass definition we can obtain the MS-masses more accurately than 
from the pole mass scheme. 

In the next section we will discuss the quantum mechanical sum rules. In 
section 3 we shall present the contributions from the threshold expansion in 
the framework of NRQCD. Here we will also define the potential-subtracted 



3 



mass. The perturbative expansion will be derived in the following section. 
In section 5 and 6 we will discuss the nonperturbative contributions and the 
phenomenological spectral function. Then we shall explain the reconstruction 
of the spectral density. In the numerical analysis we will obtain the pole- and 
MS-masses from analyses in the pole- and PS-mass scheme respectively. The 
origin of different contributions to the error will be carefully investigated. 
After a comparison to other mass determinations we shall conclude with a 
summary and an outlook. 



2 Quantum mechanical sum rules for the Cou- 
lomb potential^ 

Before studying the full quantum field theory case in detail, it will be in- 
structive to first investigate the corresponding quantum mechanical system. 
Since it is possible to describe the system analytically, one can obtain a 
clearer picture of the structure of the method and the behaviour of the dif- 
ferent contributions. Let us consider a system of two particles separated by 
a distance x. The Schrodinger equation for stationary states takes the form 



V^(x) = Eip{:>i) 



(5) 



where // represents the reduced mass of the system. The Green's function is 
constructed with help of the resolvent operator G{E) = (H — E)~^. 

By introducing a full set of intermediate states we obtain the phenomeno- 
logical side of the Green's function in position space: 



^a(x)C(y) 



En — E 



le 



+ / dE 



. p(x,y;i?') 



E' -E 



le 



(6) 



The sum runs over the discrete part of the spectrum, ipa being the cigen- 
function to the eigenvalue E^^.. The integral is taken over the continuum part 
with the spectral density p(x, y; E). By taking the derivative at x = y = 0, 
one can define a physical correlation function 



M{E) 



=y=o 



E 



{E^-E- te) 



+ 



jdE' 



P{E') 



{E'-E- ley 



(7) 



^The author would like to thank Matthias Jamin who has initiated the investigations 
on the quantum mechanical sum rules and has contributed a substantial part to the de- 
velopment of these sum rules. 
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Via a dispersion relation, the spectral density is related to the imaginary 
part of the Green's function p{E) = ImG(0, 0; E)/n. The second line of ((Tj), 
constituting the phenomenological part, can be compared to a perturbative 
expansion of M{E) thus representing the fundamental equation for the quan- 
tum mechanical sum rules [Tl) | l2() p 21j. So far the discussion has been general. 
Now we turn our attention to the Coulomb potential: 

V{r) = ~-, E^ = -^, \^^i^{0)\' = 6io6mo^ . (8) 

In order to improve the predictive power it is convenient to formulate the 
sum rules in the framework of Borel or moment sum rules. Both cases will 
give valuable information: in the Borel sum rules we can relate M[E) to the 
retarded propagator and in this way obtain explicit analytic expressions for 
the perturbative expansion. In the moment sum rules we can derive an equa- 
tion for the mass and investigate the behaviour of the pole and continuum 
contributions for different values of n and ^. 



2.1 Borel sum rules 

In the Borel sum rules one applies a Borel transformation Br to the correla- 
tion function. It is defined by 

i-E)^ 1 -E , , 
Br= m ^— V-7--T7, - = fixed. 9 

This enhances the importance of the ground state and improves the pertur- 
bative expansion. For the sum rule we will need the application of Br on 
functions of the form [x — E)~°': 

1 

7 t;^ = T^e""" . (10) 

Now we define the Borel-transformed correlator as 



M(r) = ^S.M(E) = (0|e-^^|0) 



El^"(0)l'e-''"^+ / dE' p{E')e-'^ ^ . (11) 
^ Jo 



Since the higher states are exponentially suppressed, the dominance of the 
ground state contribution for large r is clearly improved. Let us now recall 
the definition of the retarded propagator: 

ir(xf,t;;xi,t,) = (xf|e-'^(*/-*»)|xi) , (12) 
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with tf > ti and by analytic continuation one concludes immediately: 



M{t) = K{0,-iT;0,0) . (13) 

The known solution can be written as a product of the lowest order per- 
turbation theory Mo(r) and a function -^(7) of the dimensionless variable 
7(r) m 




with 



00 -| 

F(7) = l + 2v^7 + — 72 + 4v/^73^— e^'/"'(l + erf(7/n)) (14) 

n=l 

with the function erf (x) = (2/y^) /q^ exp(— t^) dt. From this formula one can 
directly deduce the perturbative expansion as a power series in small 7('r): 

/ 27r2 
M^'^\t) = Mo(r) l + 2v/^7 + — 72 + 4v/^73C(3) 

+87^C(4) + 4v/^7'C(5) + 0(7^)^ . (15) 

On the other hand, M(r) represents the sum of the poles and the continuum 
(|TT|l . By inserting the energy levels and wave functions from eq. (jHI) we can 
write the pole contributions as 

00 1 

M^°'^(r) = Mo(r)8v^7'E-Te^'^"' 

= Mo(r) (8v^7'C(3) + SV^^'d^) + 0(7')) , (16) 

where the second line shows the expansion for small 7. The continuum 
contribution can be obtained from the difference of eqs. ()15|1 and (fTT)|l . Let 
us now assume a perturbative calculation order by order in the coupling 
constant. One would successively obtain the different orders in 7". Up to 
NNLO only the continuum contribution shows up. Then, at order 7^, one 
gets a contribution from the poles ~ 8^/^T'~f^(^{3) which is partly cancelled by 
the continuum contribution ~ — 4-y7r73(^(3); likewise for higher orders. The 
poles, starting from the order 7^, can thus be expected to be suppressed in 
the numerical evaluation. However, this conclusion is premature. As can be 
seen from eq. (fTB|). the expression for M^°'^ contains an exponential in 7. 
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Whereas for very small 7 the poles can be safely neglected, for large 7 these 
contributions grow exponentially and can far exceed the continuum part. 

Let us finally note that the perturbative expansion is not only an expan- 
sion in a, but in 7 = a^J yUT/2 and thus depends on the Borel parameter 
r. The perturbative expansion and the relative size of the poles and the 
continuum will therefore depend on the sum rule window for r chosen in the 
analysis. This behaviour is better discussed in the context of moment sum 
rules which closely resemble the field theory case at hand. 



2.2 Moment sum rules 

In the quantum mechanical Coulomb problem we define the moments as 

Mn[0 = G{E) , e > 1 . (17) 



n\ n'^a \ dE 



E=(Ei 



Again, we allow for an arbitrary evaluation point ^. A natural scale is given 
by the lowest bound state energy Ei = —fia'^/2. The parameter ^ has been 
defined somewhat different as compared to eq. (jH)). The derivatives must 
be taken in an energy region below the poles where the Green's function is 
purely real. Here, we must therefore choose ^ > 1 whereas in eq. (jHl) = 
already represents a perturbative region. Solving the relevant Schrodinger 
equation, the radial Green's function is found to be 



G{r, 0; E) = ^e-'=T(l - A)f/(1 - A, 2; 2kr) 

TT 



with the variables 



A = ^ and k = ^J-2fx{E + ie) . (19) 

rC 



U{a,P;z) is the confluent hypergeometric function. G{r,0;E) is singular in 
the limit r — 0, but the moments are finite and the first moment Mi(^) is 
found to be 

Mi(0 = A (1 + 2A + 2AV(1 - A)) |^^^^^ . (20) 

The evaluation point E = C,Ei translates into A = 1/ y^. In this theoretical 
expression, all powers of A are resummed. Like the parameter 7(r) in the 
Borel sum rules, here the true parameter for the perturbative expansion is 



A = a/y—2E/fi. The higher moments can be derived from a recursion 
relation: 

A3 d 



n d\ 
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To derive the phenomenological parametrisation of the sum rules we need 
the spectral density for positive energies which can be obtained from 

p{E) = -lm G {0,0; E). (22) 

TT 

Im G{r, 0; E) is finite in the limit r — > since it is a physical quantity and 
gives the Sommerfeld factor 

p{E) = — , E>0. (23) 

Putting everything together, we obtain the phenomenological parametrisa- 
tion 



M„^''-(0 = M„^°'-(0 + °"*(0 = E 

POO 

+ dx —r- 

Jo ^ 



2 



n+1 



(24) 



Equating this quantity to the theoretical side ()20I21|) establishes the sum 
rules. Using the exact formulas, this of course represents nothing but an 
identity. The method comes into play when only limited information on 
either part is available. It can then, for instance, be used to extract the 
lowest bound state energy by solving for Ei as the higher bound states are 
strongly suppressed by the factor in eq. (j23)- Likewise, assuming that 
the energy levels are known, one can solve for the mass n since the dominant 
dependence on the mass originates from Ei: 



El 



1 

n + l ^ 



]\/[Theo _ ]\/[Cont _ j^HigherPoles j 2 



(25) 



Since the higher poles ]\^^«9^erPoies ^Qj^^ain E^. and fi as parameters, eq. (^3)) 
represents a fast converging self-consistency equation for fi. Here we shall 
not discuss the applications of this sum rule, but finally turn our attention 
to the behaviour of the pole and continuum contributions. 

To this end, we again take a closer look to eq. (j24|) . Both contributions 
depend in a similar way on n and ^. But whereas the integration runs over 
positive values of x > 0, the energy levels E^ are located below threshold. Let 
us investigate the behaviour of the moments on n for fixed ^. For low values 
of n the higher poles and high energy part of the continuum integration can 
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e = 


= 4 






^ = 


10 




n 


1 


3 


5 


7 


1 


3 


5 


7 




0.45 


0.83 


0.94 


0.97 


0.16 


0.49 


0.67 


0.78 




0.55 


0.17 


0.06 


0.03 


0.84 


0.51 


0.33 


0.22 






n = 


= 3 






n = 


7 






2 


4 


10 


50 


2 


4 


10 


50 




0.98 


0.83 


0.49 


0.10 


0.9998 


0.97 


0.78 


0.26 




0.02 


0.17 


0.51 


0.90 


0.0002 


0.03 


0.22 


0.74 



Table 1: Relative size of pole and continuum contributions to the theoretical 
moments M'I^'^° = M^°^^^ + Mf^""* for different values of n at fixed f and for 
different at fixed n. 



have a significant influence. When we proceed to larger n we enhance the 
threshold region in the continuum integration and the lowest pole in the sum. 
Taking now n fixed we see that the variation of ^ can drastically change the 
relative size of both parts: Since the moments have a singularity at ^ = 1 from 
El, values of ^ only slightly larger than one will lead to a complete dominance 
of the first bound state on the sum rules. Larger ^ enhance the higher bound 
states and also the continuum part gets more and more important. The 
results are summarised in tabled We have depicted the pole and continuum 
contributions for different values of n and S,- One should keep in mind that 
the relative size of the pole and continuum contributions does not directly 
depend on the physical system under investigation but rather on the values 
of n and ^ chosen for the analysis; namely the poles can depend strongly 
on these parameters. However, it is important to note that the main mass 
dependence originates from the first bound states since the continuum is 
largely independent of the mass. This remains true even in a region where 
the continuum dominates the moments. Therefore, to obtain good accuracy 
when extracting e.g. the ground state energy or the mass it is advantageous 
to use low ^ and high n. Then the contribution from Ei will dominate the 
sum rules. Unfortunately, in this region also the perturbative expansion 
converges more slowly. The perturbative series behaves better for lower n 
and also for higher ^ since this parameter enters directly in the expansion 
variable A = 1/a/^. Therefore, in practical applications where the exact 
solutions are not known, one must carefully choose a range of values for n and 
^ such that the theoretical calculation can be reliably trusted without loosing 
sensitivity on the parameters one would like to extract. These considerations 
will be made more explicit in the numerical analysis. Now we will discuss 
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the Coulomb contributions in the full field theory system. 



3 Coulomb resummation 

The theory of NRQCD provides a consistent framework to treat the problem 
of heavy quark-antiquark production close to threshold. The contributions 
can be described by a nonrelativistic Schrodinger equation and systematically 
calculated in time-independent perturbation theory (TIPT) ^Hl- At NNLO 
the factorised formulation has first been shown in The correlator is 

expressed in terms of a Green's function G{k) = G(0, 0, k): 

N f 4fc^ \ 

^^'^ = 2W^ + ^Gc(A:) , (26) 



where N,. is the number of colours, k = — s/4 and M represents the 

pole mass ^3]- First we will present the method in the pole mass scheme 
and afterwards discuss the PS-scheme. The constant Ch{ois) is a perturbative 
coefficient needed for the matching between the full and the nonrelativistic 
theory. It naturally depends on the hard scale and is given by ^24jj 

2 

n f \ 1 An '^sifJ'hard) . M2)r^ f OisifJ'hard)] 
'~^h[<^s) = i — 4L/ir r'-^h '^F 



TT \ TX 



Stt^ 1797r2\^ /44 Ati"^ 11 \^ 
+2iM—]+-Hlc, + C,)lJ^\, (27) 



where Cp = 4/3, Ca = 3, T = 1/2 and = 11 — In^j'i. Gc{k) represents 
the Coulomb Green's function and reads 



Gc{k) = - 



Air 



k . f k \ ^ / CpttsM 

+ In + + ^ 1 



CpasM XlJ'fac/ V 2k 



(28) 

The contributions from NRQCD are summarised in the potential. The 
Green's function obeys the corresponding Schrodinger equation 



( - ^ + Vcix) + AV{x) + ^)g(x, y, k) = 6^'\^ - y) . 



(29) 
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Here Vc{x) = —Cpag/l^l represents the Coulomb potential and AV{x) con- 
tains the NLO and NNLO corrections. The explicit form of the potential is 
given in the appendix. The full Green's function can be derived from TIPT 
and to first order one obtains 



To be consistent to NNLO one also has to apply second order perturbation 
theory to the one-loop potential. Details about this procedure can be found 
in ^3^1 • To calculate the moments from the Green's function we will 
directly perform the derivatives at s = — 4M^^ according to eq. ©. Since 
the Green's function is known analytically as a function of A; = k{s), 
this can be done numerically. In this way we take advantage of the fact 
that the perturbative expansion parameter depends on the evaluation point. 
The expansion of the moments shows the same behaviour as has already been 
discussed in the quantum mechanical sum rules ()20I2H) . There the expansion 
parameter of the moments is A = 1/ and so higher values for ^ improve 
the perturbative series. The resulting moments include both the pole and 
the continuum contributions. 

The moments depend on three scales: The hard scale fihard ~ M enters in 
the coefficient Ch of eq. (j^ . This scale is also needed for the perturbative 
expansion which will be discussed in the next section. The soft scale fisoft ~ 
Mv is a typical scale for nonrelativistic processes and the relevant scale for 
the expansion of the Green's function. Furthermore, the factorisation scale 
l^fac separates the contributions of large and small momenta and plays the 
role of an infrared cutoff. As we perform the calculation only up to NNLO 
we are left with a residual dependence on these three scales. In fact, the 
dependence of the mass on the scales, especially on fisoft, is rather large and 
will give the dominant source of the error. To obtain the central values for 
the masses we will use a set of values for fisoft, fJ'fac and fihard according to the 
physical expectations from the charmonium and bottomium systems. The 
error will then be estimated by allowing for sufficiently large variations of 
these scales. 

Though the full theoretical moments from resummation can thus be de- 
termined, we are also interested in the pole and continuum contributions 
separately. First, we want to analyse them independently and estimate their 
contribution to the error. Second, in our numerical analysis we will recon- 
struct the spectral density above threshold. At low velocities it is given by 
the imaginary part of the nonrelativistic Green's function. 

In principle, the expressions for the energies and decay widths of the poles 



G(0,0,A;) 
AG(0,0,A;) 




(30) 
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have been calculated at NNLO. One could then deduce their contribution to 
the moments as in eq. (j23)- But in this method the contributions have to 
be calculated near threshold and thus show large corrections already for the 
bottomium and cannot be trusted for the charmonium. 

Therefore we will choose a different method of evaluation f^lESl- By 
using a dispersion relation, we derive the continuum from the imaginary 
part of the correlator. From the difference we can then obtain the pole 
contributions: 



Nevertheless, for values of n and ^ used in our analysis, the poles will give 
the largest contribution to the theoretical moments and thus the dependence 
on the scales will remain relatively strong. In the numerical analysis we will 
give a detailed account on the size and behaviour of these contributions. 

Now we investigate the spectral density from the continuum part. We 
discuss the charmonium system since the differences in the expansion can be 
seen more clearly than in the bottomium which shows a faster convergence. 
In figure Q we have displayed the spectral density times the weight factor 
for the different orders. The area under the spectral density is directly pro- 
portional to the moments. We have chosen moments and scales typical for 
the numerical analysis, n = 5, ^ = 0.5, f^soft = 1.1 GeV, fifac = 1.45 GeV 
and fihard = 1.75 GeV. The dotted line represents the LO, the dashed line 
includes the NLO and the solid line represents the full NNLO result. One 
can see that the expansion of the moments converges well in the low velocity 
region. But when higher v are used, resummation is not capable to incorpo- 
rate the correct high energy behaviour. Therefore, when we reconstruct the 
complete spectral density, we will use the resummed spectral density only 
below a separation velocity v < Vgep where the expansion can be trusted. 

Part of the large corrections to the potential is not inherent to the bound 
state system but to the definition of the pole mass. In [26, ,10] it was observed 
that the long distance sensitivity in the coordinate space potential cancels 
to all orders in perturbation theory with the long distance sensitivity in the 
pole mass. Therefore a new mass definition has been proposed, the potential- 
subtracted mass mps, where the potential below a separation scale fisep is 
subtracted: 



Mi 



Poles 



n 




(31) 




(32) 
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Figure 1: Resummed spectral density times the weight factor corresponding 
to n = 5 and ^ = 0.5 in LO (dotted), NLO (dashed) and NNLO (sohd) at 
typical scales for the charmonium ^soft = 1-1 GeV, /i/ac = 1-45 GeV and 
/Uw = 1.75 GeV. 

The subtracted potential V{r, Usep) is then defined by 

V{r, fisep) = V{r) + 26m{fisep) ■ (33) 

When using these definitions the Schrodinger equation for the Green's func- 
tion takes it usual form where the pole mass is substituted by the PS-mass 
and the potential by V{r, figep) [IHl- Since the renormalon contributions have 
been subtracted from the potential, the convergence of the expansion is im- 
proved and the strength of the potential is reduced. The developed methods 
for the evaluation of the Green's function can again be employed and the 
additional contribution is absorbed in a shift of the energy. 

We then express all moments in terms of the PS-mass and perform the 
analysis for rrips. The PS-mass is perturbatively related to the MS-mass: 

ttV m J \ TT J \ m 

Wl{m, IJsep) \ ^ ( (^s{m) \^ (k -C /^'^ep W2(m,/igep) \ 

4 / v^/\ ^ m 16 y' 
Cf, ^2 = 13.443 - 1.041n/, 

190.595 - 26.655n/ + 0.653nJ, (34) 
where m = '^ms('^ms) MS-mass evaluated at its own scale. The 



"^ps(/isep) 



ki 

h 
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functions Wi and W2 can be found in the appendix where also a more complete 
list of formulas to the PS-mass is given. 

The definition of the PS-mass and its relation to the MS-mass depends 
on fisep- This scale must be taken large enough to guarantee a perturbative 
relation between the masses. At the same time it should be chosen smaller 
than a typical nonrelativistic scale as not to affect the threshold behaviour: 

AqcD < fJ'sep < M -v. (35) 

In the numerical analysis we will see that the use of the PS-mass improves 
the error of the MS-mass. To estimate the error on figep, we shall also vary 
this scale in appropriate ranges. 

4 Perturbative expansion 

The perturbative spectral function can be expanded in powers of 

the strong coupling constant a = as/n, 

rP--\s) = {s) + a (s) + R^^^ (s) + 0{a') . (36) 

From this expression the corresponding moments Ain'^^'^iO can be calculated 
via the integral of eq. For the first two terms, analytic expressions are 
available j2Z] and they read 

+2 ln(l -p))]- 4t;(3 - v^)(\n{l + p) + 2 ln(l - p)) 

-i(l - v){33 - 39v - 17v^ + 7v^) ln(p) + ^v{5 - 3v^) ,(37) 

where p = (1 — f)/(l + f) and Li2(z) is the dilogarithmic function. The 
corresponding formulas for 11'^'') and fl'-^-* can, for instance, be found in pHll29j . 
R^'^\s) is still not fully known analytically. We employ a method based 
on Pade-approximants to construct the spectral density in the full energy 
range [SHIE]- It uses available information around = 0, at threshold 
and in the high energy region. It has the advantage that it gives a good 
description until relatively close to threshold. In this region the moments 
show a strong variation for relatively small changes of the mass. A pure high 
energy expansion would only be valid for large values of the velocity and a 
matching between the threshold and the perturbative region would be less 
reliable. 
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Figure 2: Perturbative spectral density times the weight factor corresponding 
to n = 5 and ^ = 0.5 in LO (dotted), NLO (dashed) and NNLO (sohd) for 
fihard = 1.75 GeV. 



To illustrate the perturbative convergence we compare the different orders 
for the charmonium as in the last section for the resummed spectral density. 
In figure 121 we have displayed the spectral density times the weight factor, 
again for values ofn = 5, ^ = 0.5 and fihard = 1-75 GeV. The expansion 
converges well in the high velocity region. As we approach lower v, the 
expansion cannot be trusted since singular terms in v appear which have to 
be resummed. These are included in the resummed spectral density which 
sums up terms of order for n > and k = 1,2,3. The leading 

term in the perturbative spectral density at NNLO has a singular behaviour 
~ a^/f, but its contribution to the moments remains finite since the weight 
function contains a factor of v. Consequently, the graph at NNLO in figure 
121 starts with a constant set-off at f = 0. 

As will be explained in more detail in section 7, for the perturbative 
moments we will therefore mainly use the spectral density above a separation 
velocity v > Vsep with Vsep ~ 0.4. In table |21 we compare the behaviour of 
the moments for different values of C, and Vgep- The higher Vgep and ^ one 
chooses the more one approaches the perturbative region and the expansion 
improves. For typical values of the analysis, Vgep = 0.4 and ^ = 0.5, the 
convergence is under good control. 

To calculate the moments in the PS-scheme we can use the same integra- 
tion formula ^ as in the pole mass scheme, but now the spectral density 

R{s) is evaluated at the velocity v = Jl — Arripg/s and the start of the 



15 







e = 


= 











0.2 


0.4 


0.6 




LO 


0.31 


0.29 


0.18 


0.055 


^Pert 


NLO 




U.oU 


U.oz 


n no 1 




NNLO 


1 on 


U.o ( 


U.oo 


n no c: 






e = 


0.5 











0.2 


0.4 


0.6 




LO 


0.046 


0.044 


0.033 


0.014 


^Pert 


NLO 


0.096 


0.085 


0.054 


0.020 




NNLO 


0.15 


0.12 


0.064 


0.021 






i-- 


= 1 






Vsep 





0.2 


0.4 


0.6 




LO 


0.012 


0.011 


0.0092 


0.0048 


j^Pert 


NLO 


0.023 


0.021 


0.015 


0.0068 




NNLO 


0.034 


0.027 


0.017 


0.0070 








= 2 











0.2 


0.4 


0.6 




LO 


0.0017 


0.0017 


0.0015 


0.00093 


j^Pert 


NLO 


0.0030 


0.0028 


0.0023 


0.0013 




NNLO 


0.0041 


0.0035 


0.0025 


0.0013 



Table 2: Perturbative moments at LO, NLO and NNLO with jJ-hard = 1-75 
GeV and n = 5 for different values of ^. The moments are calculated only 
from the perturbative spectral density above v > Vgep- 
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Figure 3: Diagrams at 0{al) from the light quark correlator contributing 
to heavy quark production. 



integration Vgep must be transformed to this scheme as well. 

At 0(a2) it is no longer true that heavy quark production originates ex- 
clusively from the heavy quark correlator. Also the light quark correlator 
includes a four fermion cut with a heavy quark pair radiated off the light 
quarks |321- The corresponding diagrams are shown in fig. IHl When these 
contributions are included in the measurements, which of course depends 
on the experimental setup, they should be considered in the theoretical side 
as well. The diagrams have been calculated in ref. jHH]- The resulting ex- 
pressions can be split in such a way that they allow to introduce Coulomb 
resummation effects of the heavy quark pair in a straightforward way. How- 
ever, as was discussed in [32], the heavy quark pair is produced in a colour 
octet state from the gluon splitting. In this case the potential becomes re- 
pulsive and the cross section decreases close to threshold. For high energies 
the diagram gives the same contribution as the diagram with the light and 
heavy quark lines interchanged. Since the main contribution of this diagram 
comes from the perturbative region, its O(a^) contribution to the moments 
is suppressed and has a typical relative size of 2 ■ 10~^. The shift in the 
final value for the MS-bottom quark mass then amounts to Amb{mi,) fa 9 
keV and can be safely neglected within the uncertainty of this analysis. To 
conclude the theoretical side of the correlator we now discuss the condensate 
contributions. 

5 Condensate contributions 

The nonperturbative effects on the vacuum correlator are parametrised by 
the condensates. The leading term is the gluon condensate contribution. It 
has been calculated up to next-to-leading order [31] and can be written as 
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follows: 

n^r(^) = {cPAs) + aCPAs) + Oia')) . (38) 

The analytic form of the functions Cpp{s) and Cpp{s) can be found in 
The results have been calculated in dimensional regularisation with the pole 
mass M. Using a different mass prescription like the PS-mass, Cpp{s) has 
to change accordingly. In our analysis below, we employ a value of (aFF) = 
0.024 ± 0.012 GeV^ for the gluon condensate. 

Furthermore, in j^EHlEZI the dimension 6 and 8 condensate contribu- 
tions have been calculated. However, this has been done only for moments 
at = and therefore only in this case we take them into account. For 
typical values of charmonium scales, fihard = M = 1.75 GeV, their contribu- 
tion is 5%-10% of the leading gluon condensate. In fact, as will be shown in 
the numerical analysis, the absolute contribution of the condensates to the 
full theoretical moments is small, both for the upsilon and the charmonium. 
Whereas former sum rule analyses for the charmonium have emphasised the 
significance of these nonperturbative contributions, their relative suppres- 
sion in this work is due to three reasons: First, the absolute value of the 
theoretical moments increases from the Coulomb resummation. Then we 
evaluate the moments at larger ^ and smaller n where the nonperturbative 
contributions are relatively small. Finally, since we obtain a larger pole mass 
than former analyses, the condensates, starting with a power of are 
suppressed further. 



6 Phenomenological spectral function 

Experimentally, the six lowest lying ip— and T— resonances have been ob- 
served. To obtain the phenomenological content of the spectral function we 
use the narrow-width approximation for the resonances 

Rk{s) = =^ TkiA, Tfc e+e-) 6{s - El) , (39) 

where Qc,6 represents the electric charge of the charm or bottom quark, a 
denotes the running QED coupling evaluated at a scale around the reso- 
nance mass. For the charm this corresponds to the fine structure constant 
a = 1/137.04 whereas for the bottom widths the Review of Particle Proper- 
ties has used = 1.07a^ and we will do so accordingly. The narrow-width 
approximation provides an excellent description of these states since the full 
hadronic widths are much smaller than the masses. The values for the masses 
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k 


1 


2 


3 


Ek [GeV] 


3.097 


3.686 


3.770 


r, [keV] 


5.26 ±0.37 


2.12 ±0.18 


0.24 ±0.05 


k 


4 


5 


6 


Ek [GeV] 


4.040 


4.159 


4.415 


Tk [keV] 


0.75 ±0.15 


0.77 ±0.23 


0.47 ±0.10 



Table 3: Masses and electronic widths of the first six ^/'fc-resonances. 



k 


1 


2 


3 


Ek [GeV] 


9.460 


10.023 


10.355 


Tk [keV] 


1.32 ±0.07 


0.52 ±0.04 


0.48 ±0.08 


k 


4 


5 


6 


Eu [GeV] 


10.580 


10.865 


11.019 


r^. [keV] 


0.25 ±0.03 


0.31 ±0.07 


0.13 ±0.03 



Table 4: Masses and electronic widths of the first six T^-resonances. 



and electronic widths are collected in tables El and lU For our numerical anal- 
ysis the errors on the masses can be safely neglected and have thus not been 
listed. It should be kept in mind that the moments from the experimen- 
tal resonances are not identical to the ones obtained from the poles of the 
Green's function in section 3 which represent a summation of a special kind 
of theoretical contributions. 

For the upsilon system, the hadronic continuum is not measured with 
sufficient accuracy so we use the assumption of quark-hadron duality and 
integrate the theoretical spectral density above a continuum threshold Sq: 



.cstr I 



M,,n 97r ^ T,,kE,,k R^' y^, , . 

m',Y c^'Ql h (Elk + ^ 4 ' {s + 4M20"+^ • ^ 

For the parametrisation of the spectral density we use the reconstructed 
spectral density which will be discussed in the next section. The 

continuum from open B production sets in at y/s = 2Mb = 10.56 GeV just 
below the 4th resonance. In the upsilon system the resonances are relatively 
dominant. The start of the continuum threshold Sq should thus in principle 
be given by the mass of the 7th resonance. Nevertheless, when we take into 
account only the first 3 resonances and a continuum threshold Sq typically 
250 MeV above the 3th resonance and compare the result to an evaluation 
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Figure 4: Rc{s) from jSHJ. 



with all 6 resonances and Sq above the 6th resonance we miss in the latter 
a contribution of 30% from the continuum. It seems natural to assume that 
this contribution originates from open B production. To account for this 
contribution we lower the value of sq to = 11.0 GeV. To estimate the 
error we vary sq between 10.8 GeV < < 11.2 GeV. To be conservative, 
in the analysis we also check the influence on the error if we remove the 
resonances above the continuum. 

Above the charmonium threshold recent measurements have improved the 
phenomenological situation significantly [38J. 85 data points have been taken 
in the region between 2.0 GeV < ^/sq < 4.8 GeV with an average precision of 
6.6%. The continuum threshold starts at ^/s = 2Mo = 3.73 GeV. From the 
measured spectral density the light quark contributions must be subtracted. 
At this energy the light quarks can be safely assumed to be massless and the 
high energy approximation jSHl provides a good description. The resulting 
spectral density is shown in figure El At energies above the data points we 
again use the reconstructed spectral density. Since the resonances 3-6 are 
well reproduced by the data, we only add the first two resonances below 
the continuum threshold. The error from the data turns out to be small 
compared to the theoretical uncertainties. 

It is interesting to compare the measured cross section to the predictions 
from quark-hadron-duality. In average, the reconstructed spectral density 
lies above the data points. This should be no surprise as the OPE demands 
an equality of the theoretical and phenomenological moments only for the full 
correlator which also includes the pole contributions. Since the lowest poles 
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are very dominant on the phenomenological side they are compensated by a 
larger theoretical spectral density for intermediate values of s. Consequently, 
one should take care when describing the phenomenological spectral density 
by the perturbative one, in particular, the choice of the integration point 
So could depend on the values of ^ and n. A more detailed description of 
the charmonium cross section and the accuracy of quark-hadron-duality is 
presented in jlHI- Now we explain how to construct a theoretical spectral 
density for the full energy range. 

7 Reconstruction of the spectral density 

Besides the contributions from the poles of the Green's function and the con- 
densates, the theoretical part of the correlator contains the spectral density 
above threshold. Now we discuss the different parts of the spectral density. 

For high velocities the spectral density is well described by the pertur- 
bative expansion. As explained in section 4 we have used a method which 
allows a good approximation until relatively close to threshold. The re- 
summed spectral density, on the other hand, gives a good description for 
low values of v, but it fails to describe the high energy part. For these rea- 
sons, we introduce a separation velocity Vsep- Above Vsep we only use the 
perturbative spectral density. Below Vsep we essentially take the resummed 
spectral density. The perturbative expansion has singular terms in v which 
are included in the resummed spectral density, but it also contains contribu- 
tions from higher powers in v which can be isolated by subtracting the dou- 
ble counted terms and these contributions will be added to the resummed 
spectral density below Vgep- In fig- El we have displayed the charmonium 
spectral density from the different contributions as a function of v. This 
representation expands the threshold region. The dotted line represents the 
perturbative expansion at NNLO. The dashed line is the resummed spec- 
tral density and the dashed-dotted line the reconstructed spectral density. 
For the charmonium system, there exists a range of intermediate values of 
V where neither the perturbative expansion nor the resummation can be 
trusted. Indeed, it can be clearly seen that the reconstructed spectral den- 
sity shows a gap at the separation velocity. Since this gap is not physical 
but a result of the mismatch between the two energy regions we can try to 
construct a more physical spectral density which interpolates smoothly be- 
tween small and large v. We can construct this interpolating spectral density 
from i^-R*^*"™ and R^'^^^ between the two velocities Vi = 0.2 and f 2 = 0.6 with 

Jointer ^ j^Resum(^^2 _ ^;2)/(y2 _ ^2) ^ R^^^^y^ - vj)/{v^ - vj). There is UO 

explicit argument for a specific choice of except that it should give a 
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Figure 5: ^/'-system: Solid line: interpolated spectral density; daslied-dotted 
line: reconstructed spectral density; dashed line: resummed spectral density; 
dotted line: perturbative spectral density. 



smooth transition between the low and high energy region. We have chosen 
a quadratic form instead of a linear one since it suppresses better the be- 
haviour of i^^e*"™ at high V and of R^'^'^^ at low v. The moments from the 
interpolating spectral density are equal to the moments of the reconstructed 
spectral density at Vsep ~ 0.4, a typical nonrelativistic velocity. However, in 
the analysis we use the reconstructed spectral density and vary Vgep between 
0.3 and 0.5 to estimate the error. 

In figure IHl we have shown the same picture, but now for typical scales 
of the upsilon system. Here the expansion behaves much better, i^^^*"*" is 
valid up to higher and i?^*^^* to lower velocities. So we obtain an overlap in 
the intermediate region and the result is independent of variations around 



8 Numerical analysis for the bottom quark 

mass 

We now perform the analysis for the charm and bottom quark masses in the 
pole and PS-scheme. Though the method of analysis will be similar in all 
four cases we discuss every case separately as each requires a certain choice 
of parameters and an independent error analysis. Since in the upsilon system 
the theoretical expansions converge better we start with the bottom quark 



22 



8|— 
7 - 
6 - 
5 - 

3 - 
2 - 
1 



0.2 



0.4 



0.6 



0.8 



Figure 6: T-system: Solid line: reconstructed spectral density; dashed line: 
resummed spectral density; dotted line: perturbative spectral density. 



mass and devote the next section to the charm quark mass. 
8.1 Pole mass scheme 

First, one has to fix the parameters on which the sum rule depends. As a 
general rule, we will choose central values for the parameters to determine 
the masses and then vary these parameters in appropriate ranges for the 
error estimate. Let us start with the values of ^ and n. Since the bottom 
quark is relatively heavy, even for ^ = the nonrelativistic and perturbative 
expansions converge reasonably well. Nevertheless, the contributions from 
the poles of the Green's function still dominate the theoretical part. To re- 
duce their influence and to spread the theoretical contributions more equally 
among the poles, the resummed spectral density and the perturbative spec- 
tral density we must choose a higher ^. However, for ^ > 1 the moments loose 
sensitivity on the mass and the error from the input parameters increases. 
Therefore we use a central value of = 0.5 and vary ^ between < ^ < 1. 
Since the relevant scale for the evaluation point is the lowest bound state 
energy, values of ^ = 0, 0.5 or 1 already correspond to well separated eval- 
uation points. In the section on the PS-scheme we will justify this choice 
numerically as well. 

As was shown in section 2, high values of n enhance the threshold region 
and the pole contributions. To keep the theoretical expansions under control 
we restrict the moments to n < 10. From the lower side, n is limited by the 
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Figure 7: Thick solid line: central pole mass; thin solid lines: Mf, for jj^soft = 
2.0 and 3.5 GeV; dashed lines: Mb for fifac = 2.0 and 5.0 GeV; dotted lines: 
Mb for fXhard = 2.5 and 10.0 GeV. 



phenomenological uncertainty. Since for n < 4 the continuum has a large 
influence on the mass we use a range of 5 < n < 10. As central values for 
our scales we have selected: 

lisoft = 2.5 GeV , fifac = 3.5 GeV , fitard = 5.0 GeV . (41) 

We set the hard scale to the central value for the pole mass. The soft scale 
should be given by the mass times a typical velocity. We choose a value 
of /iso/t = 2.5 GeV. Though one may prefer a slightly lower value of fisoft, 
the nonrelativistic expansion gets large corrections for figoft < 2.0 GeV and 
we will use this value as the lower bound in the variation of figoft- The 
factorisation scale separates the different regions and should lie between the 
two other scales. The selected scales are required for a correct description 
of the spectral density. Since this is a physical quantity the scales must be 
chosen independent of ^ and n which merely serve as an evaluation point for 
the moments. 

As discussed in the section about the phenomenological spectral density, 
we employ a continuum threshold of ^/so = 11.0 ± 0.2 GeV. We use a sep- 
aration velocity of Vgep = 0.35 and the result is independent for a choice 
around this value. In the upsilon system the contribution from the conden- 
sates is suppressed by roughly a factor of 10~^-10~^ as compared to the pole 
contributions and can be safely neglected. 

In table El we have collected the individual moments for different val- 
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n 


5 


6 


7 


8 


9 


10 


\A Poles 
■' '■n 


0.079 


0.055 


0.038 


0.026 


0.018 


0.013 


KARcstr 
■' '■n 


0.098 


0.055 


0.032 


0.019 


0.011 


0.0067 


KA Continuum 


0.047 


0.023 


0.012 


0.0059 


0.0031 


0.0016 



Table 5: Moments for different n with the parameters ^ = 0.5, fisoft = 2.5 
GeV, /i/ac = 3.5 GeV, fihm-d = 5.0 GeV, ^/so = 11.0 GeV and Vsep = 0.35. 



ues of n. are the theoretical poles of the Green's function. Ai^'^^^^ 

contains the moments from the reconstructed spectral density and includes 
the moments from the resummed spectral density Jl4^<^^^"^ below and the 
perturbative spectral density A4^^^^ above Vgep- There is no clear distinc- 
tion between ^A/J^esum ^^^^ M^'^^^ since Vsep can be used to shift their values. 
j^Contmuum momeuts from the continuum part of the reconstructed 

spectral density above sq. For values between 5 < n < 10 the moments 
j^Rcstr g^j^g q£ ^ similar size as the pole contributions: for n = 5 they exceed 
the poles and become smaller for n > 7. The influence of the continuum 
moments is relatively strong for small n but gets more and more suppressed 
for higher n. The thick solid line in figure [U shows the central value for the 
pole mass with the scales from eq. ()4H1 . Averaging over the mass between 
5 < n < 10 we obtain 

Mh = 4.984 GeV. (42) 

The error originates mainly from the variation of the scales. For the soft scale 
we choose a range of 2.0 GeV < fisoft < 3.5 GeV. Below fisoft < 2.0 GeV 
the pole contributions show a bad perturbative behaviour and the analysis 
becomes unstable, ^soft > 3.5 GeV would use a soft scale at an energy too 
high from physical expectation. The hard scale is varied between Mf,/2 and 
2Mb, 2.5 GeV < fihard < 10.0 GeV. As before for the central values, also the 
variation of the factorisation scale should lie between the two other scales 
and we use 2.0 GeV < /i/ac < 5.0 GeV. In figure [3 we have also plotted the 
change of the mass for a variation of these scales. The error amounts to 



2.0 GeV < fisoft < 3.5 GeV 
2.0 GeV < fifac < 5.0 GeV 
2.5 GeV < /iw < 10.0 GeV 



AMfo = 95 MeV , 
AMh = 35 MeV , 
AMb = 20 MeV . (43) 



In table IHl we have listed the dependence of Ai^"'-'^'^ and A^^^**'' on the 
scales. The soft scale has a particular large influence on the poles. In table 
13 we have confronted the LO, NLO and NNLO corrections in NRQCD for 
n = 7. Instead of the A^^*^***" we directly use J\yl^'^''^'^ for this comparison. 
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f^soft 


2.0 


2.25 


2.5 


2.75 


3.0 


3.5 


J^Rcstr 


0.050 
0.034 


0.043 
0.032 


0.038 
0.032 


0.034 
0.031 


0.032 
0.031 


0.028 
0.030 


l^fac 


2.0 


3.0 


3.5 


4.0 


4.5 


5.0 


J^Rcstr 


0.043 
0.032 


0.039 
0.032 


0.038 
0.032 


0.037 
0.032 


0.036 
0.031 


0.034 
0.031 




2.5 


4.0 


5.0 


6.0 


8.0 


10.0 


J^Rcstr 


0.038 
0.035 


0.038 
0.032 


0.038 
0.032 


0.038 
0.031 


0.039 
0.031 


0.039 
0.030 



Table 6: M^"^^' and Mf^'^'' for different fisoft with fijac = 3.5 GeV and 
l^hard = 5.0 GeV, for different /x/ac with j^soft = 2.5 GeV and i^hard = 5.0 
GeV and for different /ihard with /igoft = 2.5 GeV and /ifac — 3.5 GeV. 



^soft 




2.0 


2.25 


2.5 


2.75 


3.0 


3.5 




LO 


0.024 


0.020 


0.018 


0.015 


0.014 


0.012 


J^Poles 


NLO 


0.033 


0.029 


0.027 


0.024 


0.022 


0.019 




NNLO 


0.050 


0.043 


0.038 


0.034 


0.032 


0.028 




LO 


0.022 


0.021 


0.020 


0.019 


0.019 


0.018 


Resum 


NLO 


0.014 


0.014 


0.014 


0.014 


0.014 


0.014 




NNLO 


0.016 


0.015 


0.014 


0.014 


0.014 


0.013 



Table 7: Size of the moments from the poles and the resummed spectral 
density at LO, NLO and NNLO for different values of iigoft- 
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Source 


AM5 


Variation of fisoft 

\lfiT\rit\nT\ of 11 £ 

V CaiL ICAi U 1KJ±± J. Lai f Q^d^ 

Variation of fihard 
Threshold sq 
Experimental widths 
Variation of Aqcd 
Variation of ^ 


95 MeV 
35 MeV 
20 MeV 
20 MeV 
30 MeV 
50 MeV 
35 MeV 


Total error 


125 MeV 



Table 8: Single contributions to the error of M},. 



The analysis confirms that the expansions converge better for smaller n and 
behave worse for higher n as is expected from the more sensitive testing of 
the threshold region. The corresponding shift of the mass by going from 
NLO to NNLO amounts to AMb ^ 120 MeV. 

Now we turn our attention to the other parameters. In the pole mass 
scheme, a significant uncertainty comes from Aqcd as well. To determine 
Q;s(/i) we have used a^iMz) = 0.1181±0.002 JJ and run this value down with 
the three-loop beta function. The corresponding Aqcd for three loops and 
four massless flavours is then Aqcd = 279 ± 29 MeV. Though a.^ is known 
relatively precise, the error on the mass amounts to AMf, = 50 MeV since 
the pole mass has a strong dependence on the coupling constant. The choice 
of the continuum threshold shifts the mass at low n and gives an error of 
AMb = 20 MeV. The error from the experimentally measured decay widths 
is AMb = 30 MeV. When we employ a value of ^ = the mass decreases by 
35 MeV and rises by the same amount for ^ = 1. We have summarised the 
results in table |H1 We have checked for the correlations between the errors 
from the different input parameters and have found almost no correlation 
between the errors. This holds also true for the PS-scheme and for the 
charm quark mass analysis. Thus we add the errors quadratically and our 
flnal result for the pole mass is 

Mb = 4.984 ±0.125 GeV. (44) 

Using the three-loop relation between the pole and the MS-mass which has 
been calculated recently [1111121, we obtain mbijnb) = 4.277 ± 0.116 GeV 
for the MS-mass. However, the relation between the two masses implicitly 
includes an uncertainty of 0{Aqcd)- 

Before turning to the analysis in the PS-scheme we want to investigate 
the size of the NNNLO-corrections from NRQCD. In |131 the leading-log 



27 



term of ©(a^lna^) to the energy levels has been derived in the framework 
of potential NRQCD. The results have been confirmed in 03] where also the 
corrections to the wave function at 0(q;^ In^Og) were calculated. Recently, 
the full NNNLO corrections have been computed [130^1, but only for the 
energy levels itself. The results of 03] can be used to estimate the impact on 
the Green's function and the mass. The main contributions come from the 
two lowest bound states. With a typical soft scale of fJ-soft = 2.5 GeV for the 
Green's function, the 0(af)-contributions lower the mass by approximately 
35 MeV, AMf) = —35 MeV. We have not included this mass shift in our final 
result since these contributions represent only a part of the full NNNLO 
corrections. Furthermore, the results have only been derived for |'?/'„(0)p and 
En and not for the Green's function itself. Since the expansion for the energy 
and the wave function is not very good, the mass shift may be overestimated. 
However, it could indicate the size of the NNNLO corrections. 

8.2 Potential-subtracted mass scheme 

Here the separation scale /i^ep appears as an additional parameter which 
enters in the definition of the PS-mass (jS21)- This scale should be taken large 
enough in order to guarantee a perturbative relation to the MS-mass. On 
the other hand, it should be smaller than Mv as not to affect the threshold 
behaviour. A good value is /i^ep = 2.0 GeV and we will investigate a range of 
1.0 GeV < fisep < 3.0 GeV to check the influence on the MS-mass. In figure 
IHjwe have plotted the PS-mass as a function of n. As our central values we 
obtain 

mpsfiifJ^sep = 2 GeV) = 4.561 GeV , mb{mb) = 4.241 GeV . (45) 

Relating the pole, PS- and MS-masses we make use of the recently calculated 
three-loop result for the masses [4 1 42j . Since resummation includes all orders 
in as, this choice is more appropriate than the two-loop relation. From the 
variations of the scales we obtain 

2.0 GeV < fXsoft < 3.5 GeV : Amps,b = 90 MeV , 
2.0 GeV < fXfac < 5.0 GeV : Amps,b = 35 MeV , 
2.5 GeV < fXhard < 10.0 GeV : Amps,b = 10 MeV . (46) 

In table ini we again have collected the moments for different values of n. For 
low n, Ain'^^^^ exceeds M.^"^^^ whose influence grows for larger n. In tables 
ITUl and ^2 we have shown the scale dependence and the behaviour for the 
different orders. The large coefficients which have been found in the Coulomb 
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4.8 




Figure 8: Thick solid line: central PS- mass; thin solid lines: mps,b for 
^isoft = 2.0 and 3.5 GeV; dashed Uncs: nipsfi for /ifac — 2.0 and 5.0 GeV; 
dotted lines: mps,b for A*feard — 2.5 and 10.0 GeV. 



n 


5 


6 


7 


8 


9 


10 


\A Poles 
'« 


0.043 


0.027 


0.017 


0.010 


0.0063 


0.0039 


KA Rcstr 
'■n 


0.056 


0.028 


0.014 


0.0076 


0.0040 


0.0022 


\A Continuum 
•' '■n 


0.025 


0.011 


0.0047 


0.0021 


0.00095 


0.00044 



Table 9: Moments for different n with the parameters jdgoft = 2.5 GeV, 
l^fac = 3.5 GeV, i^hard = 5.0 GeV, = 11.0 GeV and Vgep = 0.35. 



I^soft 


2.0 


2.25 


2.5 


2.75 


3.0 


3.5 


J^Poles 

J^Rcstr 


0.022 
0.016 


0.019 
0.015 


0.017 
0.014 


0.015 
0.014 


0.014 
0.014 


0.012 
0.014 




2.0 


3.0 


3.5 


4.0 


4.5 


5.0 


^Poles 
J^Rcstr 


0.019 
0.015 


0.017 
0.015 


0.017 
0.014 


0.016 
0.014 


0.016 
0.014 


0.015 
0.014 


fJ'hard 


2.5 


4.0 


5.0 


6.0 


8.0 


10.0 


J^Potes 
J^Rcstr 


0.017 
0.016 


0.017 
0.015 


0.017 
0.014 


0.017 
0.014 


0.017 
0.014 


0.017 
0.014 



Table 10: M^°^^' and Mf^'^'' for different /dsoft with fijac = 3.5 GeV and 
IJ'hard = 5.0 GeV, for different fifac with ^soft = 2.5 GeV and fihard — 5.0 
GeV and for different Hhard with i^soft = 2.5 GeV and /ifac = 3.5 GeV. 
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I^soft 




2.0 


2.25 


2.5 


2.75 


3.0 


3.5 


J^Poles 


LO 
NLO 
NNLO 


0.011 
0.015 
0.022 


0.0089 
0.013 
0.019 


0.0077 
0.012 
0.017 


0.0067 
0.011 
0.015 


0.0061 
0.0097 
0.014 


0.0051 
0.0084 
0.012 


j^^Resum 


LO 
NLO 
NNLO 


0.010 
0.0070 
0.0077 


0.0098 
0.0070 
0.0072 


0.0094 
0.0070 
0.0068 


0.0090 
0.0070 
0.0066 


0.0087 
0.0069 
0.0064 


0.0082 
0.0068 
0.0062 



Table 11: Size of the moments from the poles and the resummed spectral 
density at LO, NLO and NNLO for different values of Hsoft- 



potential at NNLO (jSHI) lead to a mass shift of Amb{mb) ^ 90 - 100 MeV 
when going from NLO to NNLO. 

When we lower the continuum threshold to ^/so = 10.8 GeV, the mass 
decreases for small n. For ^/sq = 10.7 GeV, the mass is completely stable 
for all values of 3 < n < 12. But since the analysis is too complex to draw 
conclusions on single parameters, we will not use the argument of stability 
to fix the values of the parameters. For higher values of ^/sq = 11.2 GeV the 
mass increases for small n. The error from this variation is Amps^b = 15 MeV. 
Assuming that the continuum from open B-production can be neglected, ^/so 
should lie above the highest resonance. Keeping the threshold at = 11.0 
GeV, we must remove the 6th resonance. Then the mass increases by 6 MeV. 
If we assume a very high threshold and also remove the 5th resonance the 
mass rises by 22 MeV. The error from the experimental decay widths is of a 
similar size, Ampsfi = 25 MeV. The main part comes from the width of the 
first resonance though it is not as dominant as in the charmonium system. 
In the PS-scheme, the influence of is significantly reduced. Using the 
same Aqcd as before, we obtain Amps,b = 20 MeV for the PS-mass and 
Amb{mb) = 5 MeV for the MS-mass. From this weak dependence on the 
strong coupling constant we can not use the analysis for an estimate of a^. 

Besides the contributions from resummation, this analysis includes the 
perturbative spectral density above Vsep as well. Therefore it is interesting to 
investigate the influence of the perturbative part on the analysis. In figure 
ini we have depicted the central mass with a solid line. Then we keep all 
contributions, including the poles, the resummed spectral density and the 
phenomenological part unchanged except for the perturbative spectral den- 
sity for which we use only the lowest order. The resulting mass is shown 
as a dashed line. The analysis becomes unstable for low n since here the 
perturbative contributions play an important part. Now we remove the per- 
turbative part completely and use only the resummed spectral density. The 
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4.7 




12 



Figure 9: Solid line: central PS-mass; dashed line: perturbative contribution 
only at LO; dotted line: without perturbative contribution. 



j-^sep 


1.0 


1.5 


2.0 


2.5 


3.0 


mps,b 
mb{mb) 


4.710 
4.234 


4.631 
4.237 


4.561 
4.241 


4.497 
4.245 


4.438 
4.248 



Table 12: Change of the masses for different values of /i^ep- 



dotted line signals clearly that essential information for low n is lost. Only 
for high n, dominated by the poles, the analysis becomes more stable but 
here also the expansion of the poles behaves more badly. 

In table IT^ we have varied the separation scale ^sep between 1.0 GeV < 
f^sep < 3.0 GeV. The stability of the PS-mass as a function of n remains 
almost unchanged, but the value of the PS-mass changes strongly as its 
definition depends directly on fisep- In the relation to the MS-mass this 
variation is cancelled to such an extend that the MS-mass changes only by 7 
MeV. 

Now we want to discuss the choice of ^ in more detail. Using a higher 
the theoretical expansions converge better and the dominance of the pole 
contributions is reduced. As a result, the theoretical moments are more 
equally distributed. Thus the dependence on a single contribution like the 
poles is reduced. In this way one gets a better control over the systematic 
uncertainties in the sum rules. Table ^1 shows the MS-mass for different 
^. As a measure of the uncertainty we now investigate the change of the 
mass connected with the expansion of the poles in NRQCD. The first entry 
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c 
? 


rnjjyrni,)/ [KjH vj 


Am*-"'' /fMpVl 


Am*^-* /TMpVI 


-0.25 


4.200 


71 


114 





4.215 


62 


102 


0.5 


4.241 


54 


90 


1.0 


4.262 


49 


82 


1.5 


4.278 


46 


79 


2.0 


4.292 


45 


77 



Table 13: Change of the MS-mass. when adding the difference from NNLO- 
NLO to the poles, Aml"'\ and the difference from NNLO-LO, 



Source 


AmpsM 




Variation of fJ,soft 


90 MeV 


80 MeV 


Variation of f^fac 


35 MeV 


30 MeV 


Variation of /ihard 


10 MeV 


10 MeV 


Variation of /isep 




5 MeV 


Threshold sq 


25 MeV 


25 MeV 


Experimental error 


35 MeV 


30 MeV 


Variation of Aqcd 


20 MeV 


5 MeV 


Variation of ( 


30 MeV 


25 MeV 


Total error 


112 MeV 


98 MeV 



Table 14: Single contributions to the error of mps,h and mb{mh). 



shows the central value for the mass. Now we add to the poles of the Green's 
function the difference between the NNLO and the NLO result. The increase 
of the mass is shown in the second column. In the third column we add the 
difference between the NNLO and the LO result. Here the error on the mass 
from the expansion of the poles decreases for higher ^. If we assume this 
as a conservative error estimate, for ^ = 0.5 this error is of the same order 
as the error from the variation of the soft scale. In fact, the error from the 
scales increases for higher ^ as a result of the decreasing sensibility on the 
mass. From the viewpoint of the convergence of the series, the scale variation 
tends to underestimate the error for low ^ and to overestimate the error for 
higher ^. However, for ^ = 0.5 both estimates are consistent with each other 
and we have thus chosen this value as our default. For > 1 the better 
control over the theoretical expansions is not large enough to compensate 
for the decreasing sensibility on the mass and the increasing inffuence of the 
other contributions. The error from < ^ < 1 is Amps,b — 30 MeV. Table 
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[T^ summarises the error from all contributions for mpsfi and mi){mfj). Our 
results are 

mps^b{^lsep = 2.0) = 4.561 ± 0.112 GeV , 

rrihimf,) = 4.241 ± 0.098 GeV . (47) 

This value can be compared to the MS-mass obtained from the pole mass 
scheme, mf,(mf,) = 4.277±0.116 GeV. The central value decreases by 30 MeV. 
In the PS-scheme one has better control over the systematic uncertainties, 
reflected in an improved convergence for the theoretical expansions and in 
clear perturbative mass relation. 

9 Numerical analysis for the charm quark mass 
9.1 Pole mass scheme 

The method of analysis will follow along the same lines as for the bottom 
quark mass and we will put special emphasis on the points different in both 
analyses. First we must choose the value of ^. As for the bottomium, we 
will use ^ = 0.5. At this value the pole contributions still represent the 
dominant part. In principle one would like to choose a higher value where 
the theoretical expansions converge better. However, the contribution from 
the theoretical poles varies significantly with the scales; for > 1 the mass 
depends too strongly on these variations. Thus we again use a range of 
< < 1. Again, in the PS-scheme we will investigate the perturbative 
behaviour for different ^ in more detail. Since the perturbative expansions 
converge more slowly than for the upsilon we restrict the analysis to smaller 
values of n < 7. From the lower side, we choose n > 4 since for ^ = 0.5 
the moments at n = 3 already depend significantly on the phenomenological 
part. As central values for our scales we have selected: 

f^soft = 1.2 GeV , /iy,, = 1.45 GeV, /i/,,,, = 1.75 GeV . (48) 

The hard scale corresponds to the central value of the pole mass. For the soft 
scale we would have preferred a somewhat smaller value but then the NNLO 
corrections become large. The moments for different values of n are shown in 
table El Ai^'^^^^ are the moments from the reconstructed spectral density at 
Vsep = 0.4. At this separation velocity A^^"^***^ equal the moments from the 
interpolating spectral density which was introduced in section 7. The pole 
contributions dominate the sum rule even for small n. The condensates are 
suppressed compared to the poles and have no influence on the mass. From 
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n 


3 


4 


5 


6 


7 


8 


K4 Poles 
•' '■n 


0.65 


0.48 


0.35 


0.26 


0.19 


0.14 


\yj Rcstr 


0.41 


0.21 


0.11 


0.063 


0.036 


0.021 


A/I Continuum 


0.23 


0.099 


0.046 


0.023 


0.011 


0.0058 


KA Condensates 


-0.0033 


-0.0030 


-0.0027 


-0.0023 


-0.0019 


-0.0015 



Table 15: Moments for different n with the parameters jisoit — 1-2 GeV, 
jjifac = 1-45 GeV, Hhard = 1-75 GeV and Vggp = 0.4. 



1.9 



1.8 



1.7 



1.6 




Figure 10: Thick sohd hne: central pole mass; thin solid hnes: Mc for 
fj,soft = 1-1 and 1.35 GeV; dashed lines: Mc for /ifac — 1-2 and 1.65 GeV; 
dotted hnes: Mc for Hhard = 1-4 and 2.5 GeV. 
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I^soft 


1.1 


1.15 


1.2 


1.25 


1.3 


1.35 




0.51 


0.42 


0.36 


0.31 


0.27 


0.24 


J^RcstT 


0.12 


0.11 


0.11 


0.11 


0.10 


0.10 




1.2 


1.3 


1.4 


1.5 


1.6 


1.7 




0.46 


0.42 


0.38 


0.34 


0.29 


0.25 


J^Rcstr 


0.12 


0.11 


0.11 


0.11 


0.10 


0.096 


fJ'hard 


1.4 


1.6 


1.8 


2.0 


2.2 


2.5 




0.29 


0.33 


0.36 


0.39 


0.41 


0.43 


J\^Rcstr 


0.11 


0.11 


0.11 


0.11 


0.11 


0.11 



Table 16: M^"^^' and M^^'^'' for different /j^soft with = 1.45 GeV and 
f^hard = 1-75 GeV, for different fifac with figoft = 1-2 GeV and fihard = 1-75 
GeV and for different fj^hard with fisoft = 1-2 GeV and fifac = 1-45 GeV. 



figure averaging over 4 < n < 7, we obtain 

M, = 1.754 GeV. (49) 

The error is dominated by the variation of the scales. For values of ^soft ^ 
1.1 GeV the pole contributions get large NNLO corrections and we thus 
choose 1.1 GeV < fisoft < 1-35 GeV. For the hard scale we use a range of 
1.4 GeV < fihard < 2.5 GeV and for the factorisation scale 1.2 GeV < fijac < 
1.65 GeV. Since the convergence of the nonrelativistic expansion is not very 
good for the charmonium system, the scales cannot be chosen arbitrarily far 
away from their central values. Though the analysis is stable inside the given 
intervals, the expressions tend to become unstable for scales outside of the 
chosen ranges. The error amounts to 

1.1 GeV < fisoft < 1-35 GeV : AM^ = 90 MeV , 

1.2 GeV < fifac < 1.65 GeV : AM^ = 65 MeV , 

1.4 GeV < fihard < 2.5 GeV : AM^ = 40 MeV . (50) 

Tables IT^ and IT7I show the dependence of and A^^'^'**^ on the scales 

and of }A^°^^^ and Ji4^^^'^"^ from the different orders on fisoft respectively. 

To estimate the uncertainty from as we employ Aqcd = 329 ± 29 MeV 
which is the corresponding value for 3 flavours and 3 loops. Then the mass 
shifts by 60 MeV. Since already the lowest resonances dominate the phe- 
nomenological part, the error from the measured spectral density and exper- 
imental widths is relatively small. For ^ = the mass decreases by 60 MeV 
and increases by the same amount for ^ = 1. From tabled) we then obtain 
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H'soft 




1 1 

J.. ±. 


1 15 

J. . i (J 


1 2 


1 25 


1 3 


1 35 




LO 


0.19 


0.17 


0.15 


0.13 


0.12 


0.11 


^Poles 


NLO 


0.31 


0.28 


0.25 


0.23 


0.21 


0.19 




NNLO 


0.51 


0.42 


0.36 


0.31 


0.27 


0.24 




LO 


0.097 


0.093 


0.090 


0.087 


0.085 


0.083 


j^^Resum 


NLO 


0.043 


0.046 


0.047 


0.048 


0.049 


0.050 




NNLO 


0.045 


0.040 


0.036 


0.033 


0.031 


0.029 



Table 17: Size of the moments from the poles and the rcsummed spectral 
density at LO, NLO and NNLO for different values of jisoft- 



Source 


AM, 


Variation of i^soft 


90 MeV 


Variation of ^fac 


65 MeV 


Variation of (JLhard 


40 MeV 


Experimental cross section 


5 MeV 


Experimental widths 


20 MeV 


Variation of Vsep 


10 MeV 


Variation of Aqcd 


60 MeV 


Variation of ^ 


60 MeV 


Total error 


147 MeV 



Table 18: Single contributions to the error of M,. 



the pole mass 

Mc = 1.754 ±0.147 GeV. (51) 

This corresponds to a MS-mass of mc(mc) = 1.247 ± 0.134 GeV. Again, 
there is an 0{Aqcd) uncertainty from the perturbative relation between the 
masses so we now turn to the PS-scheme to determine the MS-mass. 

9.2 Potential-subtracted mass scheme 

As in the pole scheme, we will use ^ = 0.5 within a range of < ^ < 1. 
For the separation scale we choose Usep — 1-0 ± 0.2 GeV. This represents 
a compromise value. It is still high enough for a perturbative evaluation 

and sufficiently below the hard scale. Since in the PS-schcmc the theoretical 
expansions converge better, one can employ a lower value for the soft scale 
and we will use Hsoft = 1-1 GeV. As before, we use a range of 4 < n < 7 for 
the moments. 
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Figure 11: Thick solid line: central PS-mass; thin solid lines: mps,c for 
Hsoft = 1-0 and 1.25 GeV; dashed lines: mps^c for /W/ac = 1.2 and 1.65 GeV; 
dotted lines: mps^c for IJ'hard = 1-4 and 2.5 GeV. 



n 


3 


4 


5 


6 


7 


8 


KA Poles 
'■n 


0.22 


0.11 


0.056 


0.028 


0.014 


0.0071 


KARcstr 
'■n 


0.13 


0.043 


0.016 


0.0059 


0.0022 


0.00088 


KA Continuum 


0.057 


0.016 


0.0047 


0.0015 


0.00047 


0.00015 


\4 Condensates 


-0.0016 


-0.00097 


-0.00057 


-0.00032 


-0.00018 


-0.000096 



Table 19: Moments for different n with the parameters jjtsoft = 1-1 GeV, 
Ufac = 1.45 GeV, Hhard = 1-75 GeV and Vsep = 0.4. 



In table CHI we have shown the moments for different n. The poles repre- 
sent the dominant part of the theoretical correlator. The size of the conden- 
sates is ~ 1% of the pole contributions and also in this scheme they can be 
neglected for the analysis. 

Our central values for the PS- and MS-masses are 

"^P5,c(/^sep = 1 GeV) = 1.300 GeV, mc(mc) = 1.188 GeV . (52) 

In figure ^2 we have plotted the PS-mass and the corresponding error from 
the scales: 

1.0 GeV < fXsoft < 1.25 GeV : Amps,c = 85 MeV , 
1.2 GeV < fXfac < 1.65 GeV : Amps,c = 65 MeV , 
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I^soft 


1.0 


1.05 


1.1 


1.15 


1.2 


1.25 




0.10 


0.075 


0.057 


0.045 


0.036 


0.030 


J^Rcstr 


0.020 


0.017 


0.015 


0.013 


0.012 


0.011 




1.2 


1.3 


1.4 


1.5 


1.6 


1.7 




0.083 


0.072 


0.062 


0.052 


0.043 


0.034 


J^Rcstr 


0.019 


0.017 


0.016 


0.014 


0.012 


0.010 




1.4 


1.6 


1.8 


2.0 


2.2 


2.5 




0.040 


0.051 


0.059 


0.065 


0.070 


0.076 


^Rcstr 


0.013 


0.014 


0.015 


0.016 


0.017 


0.017 



Table 20: M^^^' and M^'''*'' for different ^soft with yifac = 1-45 GeV and 
f^hard = 1-75 GeV, for different fifac with figoft = 1-1 GeV and fihard = 1-75 
GeV and for different fj^hard with fisoft = 1.1 GeV and fifac = 1-45 GeV. 







1.0 


1.05 


1.1 


1.15 


1.2 


1.25 




LO 


0.035 


0.027 


0.022 


0.018 


0.015 


0.013 




NLO 


0.064 


0.050 


0.040 


0.033 


0.028 


0.024 




NNLO 


0.10 


0.075 


0.057 


0.045 


0.036 


0.030 




LO 


0.019 


0.016 


0.014 


0.013 


0.011 


0.010 


j^y^Resum 


NLO 


0.011 


0.010 


0.0097 


0.0092 


0.0086 


0.0081 




NNLO 


0.0096 


0.0075 


0.0061 


0.0051 


0.0043 


0.0037 



Table 21: Size of the moments from the poles and the resummed spectral 
density at LO, NLO and NNLO for different values of iisoft- 



1.4 GeV < fihard < 2.5 GeV : ^mps,c = 50 MeV . 

(53) 

Table EDI and ED show the scale dependence and the behaviour for the dif- 
ferent orders. Compared to the pole scheme, the uncertainty on is much 
improved and amounts to Amp^c = 20 MeV and Amc(mc) = 10 MeV. As 
for the bottomium we now investigate the significance of the perturbative 
contribution. The solid line in figure El shows the central mass. Then we 
just change the perturbative spectral density. The dashed line shows the LO 
result and in the dotted line we neglect the perturbative contribution com- 
pletely. For low n we loose stability for the mass though the effect is not as 
pronounced as in the bottom case since for the charmonium the poles play a 
more dominant part. 

The influence of the separation scale on the masses is shown in table |221 
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■"3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' lO 

n 

Figure 12: Solid line: central PS-mass; dashed line: perturbative contribu- 
tion only at LO; dotted line: without perturbative contribution. 





0.8 


0.9 


1.0 


1.1 


1.2 


mps,c 


1.353 


1.326 


1.300 


1.277 


1.255 


rricimc) 


1.181 


1.184 


1.188 


1.192 


1.196 



Table 22: Change of the masses for different values of jjLgep- 



The definition of the PS-mass depends directly on /i^ep, but the MS-mass 
remains very stable and changes only by Amc{mf^ = 8 MeV. 

Now we turn our attention to the choice of ^. In table EHl the MS-mass 
is depicted for different ^. As we have done in table ^1 we add to the pole 
contributions the difference from the NNLO-NLO and from NNLO-LO. The 
change in the MS-mass is shown in the third and fourth column. Since the 
poles are relatively dominant, in principle we would like to use a higher 
value of ^. But from table ESI one can see that the better expansion is 
almost compensated by the decreasing sensitivity for ^ > 1. In addition, 
the errors from the other input parameters grow. Like in the bottom case, 
the error from the scales improves for lower ^. The error from the scales is 
still larger than the estimated uncertainty from the expansion, but in order 
to be conservative we will use the larger error from the scales for our error 
estimate. The variation of ^ changes the central value for the mass as well 
and for < ^ < 1 we obtain Amps^c = 20 MeV which shows a much better 
behaviour than in the pole scheme. A summary of all contributions to the 
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/;/,.(/"<.■)/ [Gf'V] 


Avil:'>/[McV] 




-0.5 


1.161 


44 


73 





1.174 


35 


64 


0.5 


1.188 


31 


59 


1.0 


1.205 


30 


56 


1.5 


1.223 


29 


54 


2.0 


1.241 


30 


54 


3.0 


1.273 


31 


54 


4.0 


1.300 


34 


56 



Table 23: Change of the MS-mass, when adding the difference from NNLO- 
NLO to the poles, Am'f \ and the difference from NNLO-LO, Am'^\ 



Source 




Amc(mc) 


Variation of Hsoft 


85 MeV 


75 MeV 


Variation of fifac 


65 MeV 


55 MeV 


Variation of fihard 


50 MeV 


40 MeV 


Variation of fisep 




10 MeV 


Experimental cross section 


5 MeV 


5 MeV 


Experimental widths 


20 MeV 


20 MeV 


Variation of Vg^p 


10 MeV 


10 MeV 


Variation of Aqcd 


20 MeV 


10 MeV 


Variation of ^ 


20 MeV 


15 MeV 


Total error 


124 MeV 


106 MeV 



Table 24: Single contributions to the error of mps,c and mc{mc). 



40 



error is presented in table |211 Finally we obtain the masses: 

mpsAf^sep = 1.0) = 1.300 ± 0.124 GeV , 

m^{m^) = 1.188 ±0.106 GeV. (54) 

This value is 59 MeV lower than the central value from the pole scheme. This 
is no surprise since in the pole scheme the theoretical contributions have large 
perturbative corrections and the relation between the masses contains large 
uncertainties as well. 

In our previous work on the charm quark mass ^Zj we obtained a slightly 
higher charm quark mass of mc{mc) = 1.23 ±0.09 GeV. In this work we have 
chosen an evaluation point of ^ = 0.5. In the theoretical QCD calculation of 
the moments in the PS-scheme we have now set the start of the continuum in 
eq. (jHT|) according to the pole mass which is the appropriate mass definition 
for free quark production. This leads to a reduction of rric by 25 MeV. On 
the phenomenological side we have included the BES data. This gives a 
better control in the region between 3.8 and 4.6 GeV where the assumptions 
of quark-hadron-duality cannot be expected to work well. Finally we have 
extended the error analysis. 

10 Comparison to other mass determinations 

Now we compare our value for the charm quark mass to other determina- 
tions. In this and the next section the MS-masses are always evaluated at 
their own scale, rric = mc{mc) and nib = 'mb{mb). The basis of the charmo- 
nium sum rules was already laid in ^3I1HI. Since then, many researchers 
have extracted the charm quark mass from the sum rules. In P the pole 
mass was determined from perturbation theory to NLO resulting in a value 
of Mc = 1.46 ± 0.07 GeV. In a second investigation ^QlllS] the analysis has 
been performed in the MS-scheme with perturbation theory to NLO. Using 
the NLO relation to the pole mass the author obtains rric = 1.26 ± 0.05 GeV 
and Mc = 1.42 ± 0.03 GeV. The author has also performed an analysis using 
resummation in LO with a value of Mc = 1.45 ±0.07 GeV. In our analysis the 
increased value of the pole mass is essentially due to large Coulomb contri- 
butions which have not been included in former analyses. As a consequence, 
the error becomes larger as well. 

The charm quark mass can also be derived from direct application of 
NRQCD to hadronic bound states. The authors of jSHlEIl have studied 
the energy level of the charmonium ground state. They conclude a MS- 
mass of nic = 1.241 ± 0.015 GeV where the error is from the variation of 
as only. In ^H] a similar analysis was performed for the pole mass with 
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the result = l.SStoH GeV. In ^ NRQCD was applied to the mass 
difference between the B and D mesons with rric = 1.21 ±0.11 GeV. Further 
improvement of these determinations may be possible in the near future. 
However, these determinations face the problem that the contributions from 
NRQCD must be directly evaluated at low energies close to threshold whereas 
in the sum rules the theoretical expansions can be evaluated in a perturbative 
region. Furthermore, nonperturbative effects may have a significant impact 
on the charmonium energy levels. 

During the last years several lattice analyses have been performed for the 
charm quark mass |SSlEllESlEniEZlEH] with rather widespread results. The 
most recent one jSH] obtains rric = 1.26 ± 0.04 ± 0.12 GeV. This calculation 
was done at a lattice spacing of a ~ 0.07 fm. Though it was done in quenched 
QCD, the authors expect a minor decrease of 5% for the unquenched case. 
Here also a discussion on the previous lattice analyses can be found. Very 
recently, two preliminary results [HHtlHIl] from different lattice groups have 
been presented with = 1.314 ± 0.05 GeV and = 1.27 ± 0.05 GeV 
(statistical error only) respectively. 

Ref. [SI] has applied pseudoscalar sum rules to heavy-light quark systems. 
Values of = 1.10 ± 0.04 GeV and rub = 4.05 ± 0.06 GeV were deduced. A 
discussion on this work can be found in jH2i| . 

In some recent works the charm quark mass has been determined from 
charmonium sum rules on perturbative grounds without Coulomb resumma- 
tion. The authors of have used moment sum rules for the charmonium. 
They use moments of n 1 — 4 and argue that resummation is not necessary 
for such low values of n. On the theoretical side H(s) is calculated up to 
0(q!^) at a scale of /z = 3 GeV. This is compared to the phenomenological 
part which includes the J/ip and ip' resonances and the data from BES |HH1 
above the D threshold. They obtain = 1.30 ±0.03 GeV. A similar analy- 
sis for the bottom yields = 4.21 ± 0.05 GeV. In (23 a contour integration 
was performed to apply the Cauchy sum rules. The integral was closed at an 
energy of 5 GeV. At this scale the high energy approximation was used to 
calculate H(s) up to O(a^). As in [22], this was compared to the two lowest 
■j/^-resonances plus the continuum data [21]. The result is rric = 1.37 ± 0.09 
GeV. Very recently, the authors of [2S] obtain rric = l-289lo;o45 GeV from a 
comparison of the perturbative spectral density to continuum data. Ref. jHE] 
has presented an update of the SVZ sum rules [2J with perturbation theory 
at NNLO. They extract the gluon condensate and a charm quark mass of 
rric = 1.275 ± 0.015 GeV. 

However, in light of the present work it seems doubtful that a reliable 
determination of the charm quark mass from the charmonium system can 
be achieved by a pure perturbative evaluation without resummation as in 
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jnSHnHinSlinni • To clarify this point, let us return to the numerical analysis. 
Choosing ^ sufficiently large and n small enough, one can easily approach 
a region where the perturbative contribution represents the dominant part. 
One could then expect that the use of perturbation theory and its relative 
small scale dependence will give a reasonable approximation. Nevertheless, 
this conclusion is misleading. 

For values of ^ = 4 and n = 4 the moment from the perturbative contri- 
bution exceeds the poles by 30% (for n < 3 the phenomenological continuum 
is very dominant and the analysis becomes unstable). Nevertheless, a vari- 
ation of the soft scale gives an error of Arrtc = 110 MeV whereas the error 
from the hard scale is Am^ = 70 MeV and this error could even be improved 
by using a high energy approximation. For even higher ^ the result is similar: 
for = 6 and n = 4 the perturbative contribution is twice as large as the pole 
contributions, but the error from fisoft gives Arric = 125 MeV and from the 
error from fihard amounts to Arric = 90 MeV. Though the pole contributions 
are relatively suppressed, the mass reacts much stronger on the remaining 
uncertainties. This behaviour can already be seen in table ESI Ever higher 
^ do not improve the accuracy any more but in same way as the theoretical 
expansion improves the sensitivity on the mass is lost. Indeed, if we set the 
pole contribution to zero in our analysis, rric would drop by approximately 
300 MeV even for very high and the analysis could not be trusted any more. 

That a description without inclusion of the theoretical poles is insuffi- 
cient can already be seen from the quantum mechanical sum rules for the 
mass (j25|) . On the phenomenological side, the main dependence on the mass 
originates already from the ffist bound states, even if the continuum part 
dominates the moments. The contributions from the poles, starting from 
0{a^), must be included in the theoretical description as well to obtain a 
reliable mass determination. As was discussed in jHEl; the charmonium sys- 
tem is not well described as a Coulomb system. In particular, the expansion 
for the higher states cannot be trusted and the effective potential may differ 
from the Coulombic one. However, it is indispensable to use the terms from 
resummation for a determination of the mass. The most important contribu- 
tion to the sum rules originates already from the ground state ()24I25|) . The 
quantum mechanical sum rules show clearly that fixed-order perturbation 
theory in a system whose ground state is governed by a Coulomb-similar 
potential leads to an instable and unreliable sum rule for the mass. 

However, it should be kept in mind that a Coulomb-dominated description 
of the charmonium system stands on less &m grounds than for the upsilon 
system whose energy is sufficiently large to allow for a reliable resummation. 
This is also reflected in the fact that the relative error of rric is almost a factor 
4 larger than the relative error of rrib. 
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QCD sum rules 


-TiUXllOI s 




lilh j [Ltc V J 


V 

V IZiZil 






KPP IH7I 






MY inn 




4.20 ±0.10 


pp dSl 


4.80 ±0.06 


4.21 ±0.11 


jp isziinn] 


4.84 ±0.08 


4.19 ±0.06 


BS HO! 




4.25 ±0.08 


H II11IZI1IZ21 




4.17 ±0.05 


This work 


4.98 ±0.125 


4.24 ±0.10 


NRQCD potential 


PY [15J 


5.015^0 Q7 


^•^'3_o.03 


Bsv l^EH 




4.19 ±0.03 






4.21 ±0.09 


Lattice QCD 


A et al. [IS! 




4.35 ±0.23 


GGMR ^ 




4.26 ±0.09 



Table 25: Some references to the bottom quark mass. 



During the last years much effort has been dedicated to the determi- 
nation of the bottom quark mass. The methods which have been employed 
were mainly based on QCD sum rules for the upsilon system, NRQCD for the 
bound states or lattice QCD. We have listed some of these results in table 1^ 
A more complete list of references can be found in p. We can directly com- 
pare our work to previous sum rule analyses. As detailed discussions about 
the advances and drawbacks of these analyses can be found in [1^70^75 j , here 
we rather want to point out some interesting differences. For the comparison 
we will use the work by Hoang j T^I7H l7^ where the most extensive analysis 
has been presented. In that moment sum rules the theoretical moments were 
directly expanded for small energies around threshold. With help of a con- 
tour integration this could be used to calculate the moments via the inverse 
Laplace transform. Since the approach has focused on the nonrelativistic 
properties of the upsilon system, for a comparison we now set the perturba- 
tive contributions in our analysis to zero. For a closer comparison we choose 
a value of = where the analysis of Hoang was performed. Two effects be- 
come important: The central value of the mass decreases by 25 MeV. Second, 
at ^ = the poles are more dominant and the influence of the perturbative 
contributions is reduced. For ^ = 0, without perturbation theory and for 
7 < n < 10 we obtain = 4.17 GeV. For n < 6 the analysis becomes 
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unstable. Three differences remain. For our central mass we have used a fac- 
torisation scale of fifac = 3.5 GeV whereas Hoang has performed a scan over 
2.5 GeV < ^fac < 10.0 GeV which roughly corresponds to a central value 
of fifac = 5.0 GeV. Second, for the parametrisation of the phenomenologi- 
cal continuum Hoang has used a continuum threshold of ^/sq = 10.56 GeV 
which corresponds to the start of BB production. As discussed in chapter 6, 
we have chosen a value of =11.0 GeV to parametrise the non-resonant 
part of the spectral density. Using /i/ac = 5.0 GeV and = 10.56 GeV 
we finally arrive at rrib = 4.14 GeV. In his work, Hoang has estimated the 
effect of a finite charm quark mass to be Ami, ^ — 30 MeV. For a massless 
charm quark he obtains mj, = 4.20 GeV which is 60 MeV higher than our 
result for similar input parameters. Our error is larger than his one by a 
factor of two. Hoang has used a x^"fit with several moments. In this way he 
gets a cancellation between theoretical contributions of the x^-function. In 
our analysis we keep the mass as a function of n which serves as additional 
check for the stability of the sum rules. Furthermore, at ^ = 0.5 the mass 
reacts stronger on a variation of the parameters. But as was discussed in 
the numerical analysis, in this way the control over the analysis is improved 
and we believe that we have thus reduced the systematic uncertainties of the 
method. 

Now we want to comment on a recent work on massless contributions to 
the heavy quark correlator jTHlIZT]. Here it was shown that at O(a^) the 
correlator contains a three-gluon massless intermediate state. Its contribu- 
tion to the moments A4n{^ = 0) contains a divergent term for s = and 
n > 4. Thus the authors have concluded that the moment sum rules can 
only be reliably evaluated at n < 4. We believe that this claim is unfounded 
as was already noticed in ^ . First, up to now the sum rule analyses contain 
perturbative contributions up to O(a^) and higher order terms from resum- 
mation. At this order the 3-gluon cut was not and should not be included in 
the analysis and all quantities are well defined. Furthermore, even if the cal- 
culation was including the full O(a^) contributions, these terms must not be 
taken into account. The terms from the three gluon cut mainly correspond 
to light quark production and the divergent parts at s = are entirely due 
to light quark production as the heavy quark pair only gives a contribution 
above s > 4M^. Since they are not included in the phenomenological part 
they also should not appear in the theoretical part and must be explicitly 
subtracted from the perturbative contributions. The remaining ambiguity 
which results from the difficulty to separate the light and heavy quark pro- 
duction in the dispersion integral above s > 4M^ is a finite effect of 0{al) 
and can be completely neglected within the uncertainty of this work. 

As discussed in j^, the problem is of more general nature. Already 
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starting at 0(a;^), it is no longer true that a specific fiavour current = 
Q'j^Q contains only heavy quark production neither that all heavy quark 
production originates from this current. As was shown in chapter 4 for the 
diagrams of fig. 01 also light-light and heavy-light correlators contain heavy 
quark production. Since only the total electromagnetic current, including a 
sum over all fiavours, is a physical observable, the single fiavour production is 
not uniquely defined. The crucial point is to set up two identical quantities: 
the phenomenological and the theoretical side should be defined is such a way 
that they contain the same contributions. The higher the order this may be 
a more and more complicated task. 

11 Conclusions 

In this work we have obtained the following values for the charm and bottom 
quark masses: 

Mc = 1.75 ±0.15 GeV, = 1.19 ± 0.11 GeV , 

Mb = 4.98 ±0.125 GeV, m,, = 4.24 ± 0.10 GeV . (55) 

As in the last section we evaluate the MS-masses at their own scale, mc,6 = 
^c,b{mc,b)- Now we summarise the key features of this analysis. 

In section 2 we have presented a rather complete setup for the quantum 
mechanical sum rules in the Coulomb potential. The correlator contains 
poles below and a continuum above threshold. The poles only start with 
a power of 0{a^), but exhibit an exponential behaviour ()l(i|l in the Borel 
sum rules or a sensitive power behaviour in the moment sum rules ()24j) . 
Therefore the relative size of the poles depends strongly on either the Borel 
parameter in the Borel sum rules or on n and ^ in the moment sum rules. The 
analysis must be performed in a certain sum rule window for r or for n and 
^ respectively to guarantee a reliable theoretical calculation and sensitivity 
to the phenomenological parameters. 

In the field theory case the expansion of the Green's function is known 
up to NNLO in the framework of NRQCD. The Green's function is directly 
evaluated at Sq = —Arn?^ and in this way we avoid to sum up the energy 
levels individually where the expansion is badly convergent. The spectral 
density can be obtained from the imaginary part of the Green's function and 
the pole contributions from the difference between the full and the continuum 
result. The resulting moments depend on three scales: Hsofti f^fac and fihard- 
In particular the dependence on figoft is relatively strong and presents the 
dominant contribution to the error. Since the pole mass contains renormalon 
ambiguities we have also performed the analysis for the PS-mass which can be 
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perturbatively related to the MS-mass. Whereas the PS-mass by definition 
depends on the separation scale which was used to subtract the long-distance 
potential, this dependence cancels in the transition to the MS-mass to a large 
extend. We have then included the perturbative contributions up to O(a^). 
They are necessary to construct the spectral density above threshold for the 
full energy range and to guarantee the stability of the mass in a region of 
small n. 

One of the great virtues of the method of QCD sum rules is the analytic 
dependence on the theoretical and phcnomenological parameters. Thus we 
have investigated their importance and influence on the analysis. For the 
determination of the masses we have used central values for all parameters. 
These values were not motivated by any optimisation or stability require- 
ment, but only grounded on general considerations. Each was varied in a 
suitably large chosen window for the error estimate. Only the lower value 
of fisoft was also limited by the convergence of the nonrelativistic expansion. 
Finally, all errors have been added quadratically. 

We have set up the sum rules for an arbitrary evaluation point ^. With 
this parameter it is possible to shift the moments into a more perturbative 
region for higher ^ or into a region more sensitive to the bound state energies 
and the mass for lower ^. We have used ^ = 0.5 both for the charmonium 
and bottomium. In fact, moving away from the threshold region and loosing 
sensitivity on the mass, the scale dependence is even a bit larger than at 
^ = 0. But the convergence of the theoretical expansions is improved and 
the theoretical contributions more equally distributed among the different 
terms, in particular, the pole contributions do not play such a dominant 
role. Thus we believe that we have reduced the systematic uncertainties in 
this sum rules which cannot be accounted for by a variations of the scales. 

We would like to emphasise a remarkable property of this analysis: Once 
a particular set of (central) values for (Xsoft, l^fac and fihard is chosen, the MS- 
masses remain very stable over a large range of values for n, ^ or ^sep- In gen- 
eral, the MS-bottom quark mass changes only by Ami, = ±30 MeV from its 
central value for 4 < n < 15 and the charm quark mass by Arric — ±15 MeV 
for 3 < n < 15. Varying ^ between —0.25 < ^ < 2 the bottom quark mass 
changes by Arrih MeV and the charm quark mass by Atoc =~^lo MeV 

for —0.5 < < 2. The variation of Usep changes the PS-mass signiflcantly 
since its definition depends on Usep- But in the relation to the MS-mass this 
change is almost completely cancelled and the bottom quark mass changes 
by Arrib — ±7 MeV for 1.0 GeV < Usep < 3.0 GeV and the charm quark 
mass by Arric — ±8 MeV for 0.8 GeV < Usep < 1-2 GeV. These results are 
astonishing since the variation of these parameters corresponds to largely dif- 
ferent relative influence among the theoretical contributions. Thus we hope 



47 



that we were able to set up a consistent framework in which the physics 
of the relevant energy region, apart from the remaining uncertainties in the 
nonrelativistic and perturbative expansions, has been correctly described. 

Let us finally summarise the achieved status. In our analysis, several 
contributions seem to be under good control: The perturbative expansion, 
as it has been incorporated in our analysis, converges reasonably well. The 
condensates only give a negligible contribution to the upsilon or charmonium. 
On the phenomenological side, for the upsilon the first six resonances have 
been measured. For the non-resonant continuum part, quark-hadron duality 
has been used. In the charmonium system, the experimental situation has 
improved recently. Besides the first six resonances also the cross section be- 
tween 2.0 GeV and 4.8 GeV has been measured. Above this energy we again 
use quark-hadron duality. Nevertheless, the most important contribution is 
already given by the first two poles. 

Of decisive importance for the determination of the masses is the thresh- 
old behaviour. The method of QCD sum rules is a very powerful tool to 
extract the masses since - by the choice of n and ^ - it can react very sen- 
sitive to the threshold. Thus, large theoretical uncertainties only lead to a 
relatively small shift in the masses. The main uncertainties indeed come from 
the threshold expansion of NRQCD. The largest potential for an improve- 
ment of the analysis lies in a further understanding of this energy region. In 
particular, the knowledge of the Green's function at NNNLO could help to 
reduce the error. 

The method of QCD sum rules is based on the assumptions of quark- 
hadron duality. With the development of NRQCD is has become clear that 
the pole contributions must be included in the theoretical description for a 
correct comparison between the theoretical and phenomenological part. The 
theoretical description is based on the operator product expansion and can 
be performed in a perturbative region where all expansions converge well. 
However, this cannot be used for an ever increasing precise determination 
of the mass. Since on the phenomenological side the dependence on the 
mass originates mainly from the first bound states, the sensitivity on the 
mass decreases in a similar way as the theoretical expansions improve. In 
our analysis we have tried to balance these contributions by an appropriate 
choice of ^ and n. Without significant progress in the theoretical description 
it seems that further substantial improvement will be difficult to achieve. 
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A Potential 

In NRQCD the Green's function obeys the Schrodinger equation: 

+1^naV(x) + i.BFV(x) + ^ j G(x, y, 4) = ^l" (x - y) . (56) 

Vc{x) is the CoTilomb potential, ^NAy{x) the nonabelian part of the quark- 
antiquark potential, ^BFVix) the Breit-Fermi potential and the terms l^iV{x) 
and ^2V{x) contain the first and second order perturbative correction to the 
Coulomb potential. The exphcit expressions read: 

yc{x) = , a; = X , 

X 

^naV[x) = -CaCf 



s 



2Mx^ ' 

:(3)(x)- 

Cpas fS^ 3(Sx)2 47r 



^bfV(x) = -^*<)(x)- — (^p +-x(xp)pj+^^SL 



2Mnx= x= 3(2S^-3)^'"(x)j, 

Ai\/(a;) = Vc{x) (ai + 26o7i? + 26o ln(x/i)) , 

^2V{x) = \/c(a:)(a2 + 6^7rV3 + 47|)+ 27^(61 + 26oai) 

+ (261 + 46oai + Sbl^E) In(xAt) + Ahl W {x^x)) , (57) 

with the constants 
2nf 



bo = 11- 



3 



34 20 

bi = —Cl - —CaTuj - ACpTuf = 102 - 12.67n/ , 

o o 

31 20 
oi = ~g^^ — g""^^-^ ~ 10.33 — l.lluf , 



49 



/4343 , . TT^ 22,, A ^„ /1798 56,,,,^^ 
«2 = (^7^ + 4vr^-X + y^(3)jcl-(^ + -C(3))C^rn; 

- (f - 16C(3)) C.Tnj + (^Tn, 

= 456.75 - 66.35n/ + 1.23nJ , (58) 

where Ci? = 4/3, Ca = 3 and T = 1/2. The coefficient 02 was ffist calculated 
in [ZHl and later corrected to the above given value W^ . 



B Potential-Subtracted mass 

The PS-mass is defined by [JH] 

q|<A'sep 

mpsifJ^sep) = M - 5m{fXsep) ■ (59) 

The subtracted potential V{r,fisep) is related to the potential in momentum 
space V{q) and 6m{fisep)- 

V(r,^„p) = V-(r) + 2<Sm(^„,) , (60) 



with ai and 02 as in eq. ()58|) . Performing the Fourier transformation one 
obtains the relation between PS- and pole mass: 

mpsilJ^sep) = M 1 + a(/i)r25(/isep) + a?{^j)r%{^i, Usep) + a^{^j)rfl{^, fi^ep) 

(1)/ ^ _ ^ A^sep ^(2)/,, N _ fJ'sep WiifX, fXsep) 

^PSKH'sep) — ■> ^PSKH'^ H'sep) — —'^F^^ ^ , 

rgk/^,/^.e,) = -Cp^^^^^il^, (61) 
with a = as/iT and the functions 

WiilJ', i^sep) = ai-bo (in ^ - 2^ , 

W2{fi, fisep) = a,- {2a,bo + h) (\n ^ - 2^ + bl (in' ^ - 4 In ^ + 8^ , 



(62) 
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with bo and bi from eq. Using the three-loop result between pole and 

MS-mass one can relate the PS- and MS-mass: 



mpsifJ^sep) = rn 



1 + a(m) ( h - Cf^) + a\m) ( h - 

V m J \ m 



Wi[m, fisep) 

4 



+a (m) ks - Cf — 

y m lb J 

= Cf, fca = 13.443 - 1.041n/ , 
fcg = 190.595 - 26.655n/ + 0.653nJ , (63) 

where m = ?^ms(^ms) MS-mass evaluated at its own scale. 
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